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Abstract 

A consistent quantum mechanical calculation of partial cross sections leading to different final 
states of antiprotonic helium atom was performed. For the four-body scattering wave function, 
corresponding to the initial state, as well as for the antiprotonic helium wave function, appearing 
in the final state, adiabatic approximations were used. Further, symmetric and non-symmetric 
effective charge (SEC, NEC) approximations were introduced for the two-electron wave functions 
in the field of the two fixed charges of the He nucleus and the antiproton. Calculations were carried 
out for a wide range of antiprotonic helium states and incident energies of the antiproton below 
the first ionization threshold of the He atom. The origin of the rich low-energy structure of certain 
cross sections is discussed in detail. 
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I. INTRODUCTION 

One of the most impressive success stories of the last decade in few-body physics is the 
high precision experimental and theoretical studies of long lived states in antiprotonic helium 
(for an overview see While the energy levels have been both measured and calculated 
to an extreme precision, allowing even for improvement of numerical values of fundamental 
physical constants, some other relevant properties of these states were studied with con- 
siderably less accuracy. Among these is the formation probability of different metastable 
states, characterized by total angular momentum J and "vibrational" quantum number v, 
in the capture reaction 

p + ^He—^{^He+p)j, + e-. (1) 

The existing calculations of the capture rates of slow antiprotons in He P, 0, Oj are based 
on classical or semiclassical approaches and they mainly address the reproduction of the 
overall fraction (3%) of delayed annihilation events. Recent experimental results from the 
ASACUSA project 5|, however, yield some information on individual populations of different 
metastable states, and our aim is to perform a fully quantum mechanical calculation of the 
formation probability of different states in the capture reaction. 

The exact solution of the quantum mechanical four-body problem, underlying the re- 
action (0) is far beyond the scope of this work, and probably also of presently available 
calculational possibilities. Still, we want to make a full quantum mechanical, though ap- 
proximate, calculation of the above process. Full is meant in the sense that all degrees of 
freedom are taken explicitly into account, all the wave functions we use, are true four-body 
states. 



II. CALCULATION METHOD 



The partial cross section, leading to a specified final state (J, v) of the antiprotonic helium 
can be written as 



(2) 
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where is the exact 4-body scattering wave function corresponding to the initial state 

$5,,,K,(ri,r2,R) = $He(ri,r2)^^^e'^>^, (3) 
while the final state j^^. is taken in the form 

•^i.K/Pi, P2, R) = '^JviPi, R) e^^^^P^ . (4) 

Here are the vectors pointing from the helium nucleus to the i-ih electron, R is the vector 
between He and p, and are the Jacobian coordinates of the electrons, measured from the 
He — p center of mass: 

ri = pi + aK; a = , (5) 

rup + niHe 



while /ij and /i/ are the reduced masses in initial and final channels, respectively. In Eq. (jH)) 
$He(i"i,i"2) denotes the He ground state wave function, while in Eq. (jl)) $jt,(oi,R) is the 
antiprotonic helium final state, for which we used a Born-Oppenheimer form [alD]: 

$,.(p, R) = Yjm{R) vt''\p; R) (6) 

where '^i^' ^ (p; R) is a two-center wave function, describing the electron (ground state) 
motion in the field of two charges (Zi, Z2), separated by a fixed distance R: 

(-^Ar + f + ^) ^1?'^^V; R) = et-'^\R) ^S?'^^^; R) (7) 

while Xjv{R) is the heavy-particle relative motion wave function, corresponding to {^Hep e~) 
angular momentum J and "vibrational" quantum number v. 



(8) 



/i being the He — p reduced mass. 

The transition potential in Eq. Q is obviously the interaction of the emitted electron 
(#2) with the rest of the system: 

r2 Irs - R| Iri - ral 



^ Since in the following we shall deal only with la two-center ground states, the index la will be omitted 
throughout the paper. 
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The electron anti-symmetrization is accounted for by taking an ri r2 symmetric initial 
state wave function {S = 0) and the factor 2 in front of the cross-section reflecting the 
indistinguishability of emitted particles sl- 



The general expression Q for the cross-section, leading to a specific state (J, v) can be 
rewritten in terms of matrix elements between angular momentum eigenstates as 

aj, = 2 (27r)^ ^ /i, //^ ^^(27^ + 1) \Mjf (10) 
* Jt,i 

with 

= ( R) <jyK,Ap2)YA I Vf I ^SX*(Pi, P2, R) ), (11) 

where [ stands for vector coupling, ^^hI'^I exact scattering wave function with 

total angular momentum J^, corresponding to the initial state 

and denotes free states with definite angular momentum 

It can be seen from Eqs. fllOlllj] . that a given antiprotonic helium final state (J, t') can be 
formed from different total angular momentum states, depending on the orbital momentum 
/, carried away by the emitted electron. 

The simplest way of approximate evaluation of Eq. (j21) or (fTIH) is to use Born approxi- 
mation, replacing the exact scattering wave function '^^^e'K ^"7 its asymptotic form ^^gK- 
from Eq. Q. In order to get an idea of the feasibility of such a "full" (including all degrees 
of freedom) calculation we evaluated the cross-sections aj^ in Born approximation in a wide 
range of quantum numbers {J,v). For the He ground state wave function in this case we 
used the simplest variational form 

$//e(ri,r2) =Nexp(-(T(ri + r2)) (12) 



with cr = 27/16 taken from book j9|. In spite of the known poor quality of the Born 
approximation for slow collisions, due to the realistic final state wave function, we hoped to 



get some information at least about the relative population probabilities of different final 
states. These expectations were not confirmed, the Born cross-sections turned to be orders 
of magnitude away from the more realistic ones. The detailed results of the Born calculation 



can be found in 



id 



There are two basic drawbacks of the Born approximation for slow collisions and long- 
range forces: 

— the antiproton "feels" the interaction from the He atom, it approaches, therefore, its 
wave function in the form of a plane wave has to be modified, 

— the He electrons also "feel" the approaching antiproton, the polarization of their wave 
functions has to be taken into account. 

To meet these requirements we used an adiabatic, Born-Oppenheimer type approximation 
for the wave function 

ne,K,^'^'He(ri,r2;R)xK,(R), (13) 

where $//e(ri, R) is the ground state wave function of the He atom in the presence of a 
negative unit charge (the antiproton) at a distance R from the He nucleus: 

■HHe,p{R)^He{ri,V2]K) = e{R) ^ He{vi, ^2]^) , (14) 

112 2 1 1 1 
HHeAR) = -o^ri - T^Ara - - - - + -7 + j- w] + 



2 2 ri r2 — 1"2| ki — R| |r2 — R| 

and (R) is the antiproton scattering wave function in the adiabatic He — p potential: 

VHe-p{R) = -l+e{R). (15) 



+ VHe-p{R?j XK.(R) = ^ XK.(R). (16) 



-2 

2/i 

In this approach the most difficult task is the solution of Eq. ()14j) . the determination of 
the wave function of two interacting electrons in the field of two fixed charges. Instead of 
performing a cumbersome variational calculation, as e.g. in jlll . 12 1, we follow an approxi- 



mation scheme proposed by Briggs, Greenland, and Solov'ev (BGS) according to which 
the solution of Eq. ()14|) can be sought in the form of two single-electron two-center wave 
functions: 

$He(ri,r2;R) ^ ^(^"'^^^^(ri; R) (^(^^i'^^^)(r2; R) (17) 
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with 

(A^r + ^ + -^^] ^(^-^-)(ri; R) = e^^'^'^'^\R) ^(^-^-^(ri; R) (18) 
\ 2 r |r — R| / 

and the e{R) of Eqs. (|14l15jl is 

e{R) = + e^^^''^^^\R). (19) 

In this construction the effect of the electron-electron interaction |ri — r2|~^ in Eq. (fT^ is 
approximated by suitable choice of the effective charges {Zn, Z12, Z21, 222)- BGS suggest to 
fix the effective charges "seen" by the first electron, Zn and Z12, at the real charges of He 
and p, while those for the second one, Z21 and Z22, can be obtained from the requirement, 
that in the two limiting cases i? ^ and R 00, the ground state energies of H~ ion and 
He atom should be reproduced: 

eiR^O) = Egs{H-), e{R^oo)=Eg,iHe). (20) 

The conditions ()2n|l are fulfilled for 

Zu = 2.0, Z12 = -1.0, 

Z21 = 1.3444, Z22 = -1.1095. (21) 

For intermediate R-s b{R) is given by ()19|1 . 

As for He wave function, the two electrons in this case are treated in a non-symmetric 
way, and the wave function has to be symmetrized explicitly: 

$^^e(ri, T2; R) = N{R) [(/.(^"■^-)(ri; R) ¥;(^-'^-)(r2; R) + 



There is, however, a more symmetric realization of the BGS idea: starting with the plausible 
requirement, that the two electrons should "see" identical effective charges: Zu = Z21, 
Z22 = Z12 we still can impose the conditions ()20p for R —* and R 00, only in this case 
the e{R) will be the sum of two equal single-particle energies: 

e{R) = 2e^^^^'^^^\R). 

For this case we get 

Zu = 1.704, Z22 = -0.9776. (23) 




1 1 1 1 1 1 1 1 1 1 1 1 1 

1 2 3 4 5 6 

R, a.u. 

FIG. 1: Electronic energies e(i?) for the NEC, SEC and variational cases. 

The e{R) in this case is very similar to the previous one, maybe a httle closer to the "quasi- 
exact" variational curve. In this second case — for brevity let us call it SEC (Symmetric 
Effective Charge), in contrast to the NEC (Non-symmetric Effective Charge) case — the 
wave function is simply 

$He(ri,r2;R) = ^.(^"'^^^^(ri; R) (/.(^"■^^^^(rs; R) . (24) 

The differences between electronic en ergi es e{R) for the NEC, SEC and variational 
calculations (performed by Ahlrichs et al and Gibbs jl^) are shown on Fig. Q It is seen 
that both cases reproduces the variational results remarkably well, while SEC is practically 
indistinguishable from the more recent of them [l^ . 



For both choices ()22|1 and ()24|1 the definite total angular momentum wave function cor- 
responding to (fT^ can be written as 

*S,(ri>r2,R) = $He(ri,r2;R)^|^y;,M.(i?) , (25) 

since the Icr ground state two-center functions (y9(r; R) do not carry any total angular mo- 
mentum: they are eigenfunctions of = {Ir + L-r)^ with zero eigenvalue, even if they are 
not eigenfunctions of and separately. The function Xxii-^) satisfies the equation 

xU^) = (26) 
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NEC potentials, J, = 30..45 




R, a.u. R, a.u. 

FIG. 2: Effective He — p potentials Vj^(i?) for different Jt values. 



with the effective He — p potential 



(27) 



To solve Eq. ()2()j] numerically, first, the asymptotic form of XK^iR) has to be clarified. The 
asymptotic behavior of the Icr two-center energies can be written as 



e 



{Zl,Z2) 



(R) 



R^oo 



Z2 

R 



0{R 



(28) 



and thus using Eqs. (fT^ and (fT^ we get 



VHe-p{.R) 



^11 



■'21 



2 + Z12 + Z. 



^22 



+ 0{R-^) 



(29) 



2 R 

Dropping the irrelevant constant term from (j^^ we see, that asymptotically it corresponds 
to a Coulomb-interaction with effective charge Zas = — (2 + Z12 + ^22)- From the actual 
values of Z12 and Z22 (j2^ and (j^Hj) we can conclude, that NEC corresponds to a weak 
repulsion, while SEC — to an even weaker attraction. In reality, of course, there is no 1/R 
term in the asymptotic He — p interaction, since the He atom is neutral. 

It has to be noted, that in spite of the closeness of the NEC and SEC electron energies 
on Fig. ^ when we include the centrifugal term, the depth of the minima and the height 
of the potential barriers differ significantly (see Fig. |2I) and this fact strongly influences the 
low energy capture cross sections. 

According to ()29|1 Eq. (j26|) has to be solved with the asymptotic condition 



x'l{R) cos5jXK,)FjXr,,K,R) + s\n5jXK,)GjXr,,K,R) 



(30) 



8 



where Fj^ and Gj^ are the regular and irregular Coulomb wave functions, with Sommerfeld- 
parameter 



and 6j^{Ki) is the phase shift caused by the non-coulombic part of the potential. After 
the numerical solution of Eq. (j^Bj) with boundary conditions ()3U|) the matrix elements 
entering the formula (jlUj) for the cross sections can be calculated by numerical integration. 



III. RESULTS AND DISCUSSION 

We start the discussion of our results by the remark, that the expression for the 
matrix element Mj] in our adiabatic approximation can be rewritten as 

Mj} ~ J XJviR)SiiR; Kj)x-iiR)dR , (32) 

where xjv{R) and Xk,{R) are the He — p relative motion wave functions, introduced in 
Eqs. dH) and (j^ . respectively, while Si{R;Kf) contains "all the rest": the three poten- 
tials © integrated over electron wave functions and coordinates, angular variables of R 
and summations over intermediate quantum numbers. This representation is useful, since it 
turns out, that the basic dependence of the matrix elements on the quantum numbers and 
incident energy is contained in the two x functions, while Si{R;Kf) weakly and smoothly 
depends on its arguments with a significant decrease with increasing / — the orbital momen- 
tum of the emitted electron. For a few selected cases the three functions in the integrand of 
Eq. (j32|l are shown in Fig. El This feature of Si{R; Kf) allows the interpretation of Eq. p2| 
as a matrix element of antiproton transition from the initial state Xk^{R) into a final state 
Xjv{R) under the action of the effective potential Si{R; Kf). 

We have calculated the capture cross sections leading to different final states for an- 
tiproton energies below the first ionization threshold Eiab = 30.8 eV. The overall energy 
dependence of the NEC and SEC cross sections crjy{E) is shown in Figs. |3J El for a few 
quantum numbers from the region of expected largest capture probability. All cross-sections 
are measured in units of Og, oq being the atomic length unit. The main features of these 
curves can be summarized as follows. 

Obviously, final states with energy below the He atom ground state energy (-2.9036 a.u.) 
have a positive Q value, so they can be reached for arbitrary low antiproton energy. For 
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1 2 3 1 2 3 

R, a.u. R, au. 



FIG. 3: Functions XJv-, Xk^ and Si for different sets of parameters. 

example, on Figs. EJ |31 there are such states: with (J = 34, f = 0,1, 2), (J = 35,v = 0,1), 
and (J = 36, f = 0). States with higher energy can be excited only above their threshold 
energies; the steep rise of the cross sections above these thresholds can be clearly seen. 

Another remarkable feature of certain cross sections is their rich low energy structure. 
This is due to the repulsive barriers of the effective potentials Vj^ for < 38 — 39, as seen 
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FIG. 4: Energy dependence of the cross-sections for some {J,v) states, NEC case. 

on Fig. 121 These barriers, in general, strongly suppress the penetration of x^. (R) into the 
interior region, thus reducing the sub-barrier cross sections. For certain sub-barrier energies, 
however, there are quasi-stationary states in these potentials, when the interior wave function 
has a large amplitude, leading to sharp resonances in the cross sections. In order to clarify the 
origin of these peaks, we looked at the energy dependence of the phase shifts of Eq. (jSHl)- 
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0,20-, 
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6 8 10 12 14 8 9 10 11 12 13 14 15 

FIG. 5: Energy dependence of the cross-sections for some (J, v) states, SEC case. 

On Fig. ini we plotted the quantity d6j^{E)/dE — the so called time delay — which for 
isolated resonances is very similar to the more familiar Breit-Wigner cross section curve. It 
can be seen, that for all angular momenta Jt for which the potential has a barrier, there is a 
narrow resonance which is correlated with a corresponding peak in the capture cross section. 
A given cross section curve may contain several of these peaks, corresponding to different 
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FIG. 6: Time delays for different values of total angular momentum Jt, NEC and SEC cases. 



Jt and / values contributing to formation of a given final state, according to tlie sum in 
Eq. (fTUI) . In general, it is interesting to note, that in contrast to a common belief, the sum 
of Eq. |TO|l is not dominated by the s-electron emission (J = Jt,l = 0) term, the j9-electrons 
practically always, while the ci-electrons in certain cases contribute significantly. The reason 
for this may be, that the decrease of Si{R; Kf) with growing / could be "compensated" by 
the possibility of lower Jf values, for which the effective potentials contain less repulsion, 
thus allowing more penetration of Xr, into the interior region. 

The d6j^{E)/dE plots apart from the narrow peaks corresponding to quasi-stationary 
states, show another, much broader peak, in some cases superposed on the narrow one. 
This one is connected with specific behavior of elastic scattering when the energy is close to 
the potential maximum; it is called "orbiting" Q|. 

The behavior of the incident antiproton wave function (R) for different energy-regimes 
with respect to the barrier maximum are illustrated in Fig. [7| 

Final states with higher J, for which the relevant effective potentials have no barrier 
show a simple energy dependence: a steep rise above the threshold and then an exponential 
decay for higher energies. The exponential fall of the cross sections for increasing energies is 
characteristic for both barrier-posessing and barrier-less potentials and is due to increasingly 
rapid oscillations of Xk (R) iii the interior region which reduce the value of the integral in 

Eq. dsa). 

The quantum number dependence of certain cross sections is shown in Fig. |S1 for some 
above-barrier energies, where such a comparison makes sense. 
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0,0- 



01,0- 



R, au. R, a.u. 

FIG. 7: Incident antiproton wave functions (in arbitrary units) with their potentials. The incident 
energies (in CM, the baselines for the wave function plots) are chosen to represent different cases 
with respect to barrier maxima. 

a) - sub-barrier non-resonant energy; b) - sub-barrier resonant energy; c) - "orbiting": energy 
close to the potential maximum; d) - above-barrier energy. 

IV. CONCLUSIONS 



To our knowledge, this is the first fully quantum mechanical calculation of the process (P), 
with all degrees of freedom taken explicitly into account. The adiabatic wave functions 
used both for initial and final states seem to be reasonably realistic. The results show, 
that quantum mechanical treatment is really necessary, especially in the low-energy region, 
where barrier penetration and resonance effects are essential. The energy dependence of the 
calculated cross sections show, that the different final states {J,v) are excited with a large 
probability in a relatively narrow window of the incident antiproton energy. In principle, 
this property could be used for selective excitation of certain states. On the other hand, the 
strong energy dependence of the cross sections prevents us from making statements about 
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FIG. 8: Cross-sections for the lowest few vibrational quantum numbers v and different incident 
antiproton energies (NEC and SEC calculations) plotted against the principal quantum number 



the experimentally observable population numbers of different states since the initial energy 
distribution of the antiprotons before the capture is unknown. Even if this distribution was 
known, the observed and calculated population numbers could deviate due to collisional (or 
other) de-excitation of states in the time interval between the capture and the measurement. 
Nevertheless, we plan to make calculation of primary populations taking some trial energy 
distributions for the antiprotons. 

In the discussion of our results we deliberately did not take a stand concerning the 
NEC and SEC approximations. In general, the structure of both cross sections (energy- 
and quantum number dependence) is similar, however, SEC gives considerably larger cross 
sections, probably due to the somewhat larger attraction of the SEC effective potentials. 
We personally think, that SEC is physically more realistic, and the coincidence of SEC's 
electronic energies with those of recent variational calculation can be seen as some 
confirmation for this point of view. 
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